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Abstract. We investigate numerically the response of an atomic Bose-Einstein 
condensate to a weakly-elliptical rotating trap over a large range of rotation 
frequencies. We analyse the quadrupolar shape oscillation excited by rotation, 
and discriminate between its stable and unstable regimes. In the latter case, 
where a vortex lattice forms, we compare with experimental observations and 
find good agreement. By examining the role of thermal atoms in the process, we 
infer that the process is temperature-independent, and show how terminating the 
rotation gives control over the number of vortices in the lattice. We also study the 
case of critical rotation at the trap frequency, and observe large centre-of-mass 
oscillations of the condensate. 



PACS numbers: 03.75.Lm,03.75.Kk 
1. Introduction 

Recent experiments on dilute atomic Bose-Einstein condensates (BECs) have 
generated vortex lattices by mechanically rotating the system [fl |2j in analogy 
to the classic rotating bucket experiment in liquid Helium 0. The smooth magnetic 
trap confining the condensate, with trap frequency ou r , is deformed anisotropically 
and rotated at frequency Q. Shape oscillations of the condensate are excited, and 
the quadrupole mode plays a pivotal role in the evolution of the system: for low 
rotation frequencies, Q < 0.7cu r , a long-lived quadrupole mode has been observed 0, 
while for higher rotation frequencies Q ~ 0.7o; r , vortex lattice formation has been seen 
m0 EH- In the latter case, vortex formation is seeded by a dynamical instability of the 
quadrupole mode jSJ [7[ |SJ . Experimental results suggest that the crystallisation of the 
vortices is temperature-independent jHj. Indeed, the full process of lattice formation 
has been successfully simulated using the zero temperature Gross-Pitaevskii equation 
[Tol ITTl IT2] (although other studies found it necessary to include damping effects 
These results have also indicated that lattice formation is a two-dimensional 
effect. An experimental aim in this area is to examine the regime of fastly rotating 
BECs, and probe exotic phenomena such as vortex melting and the bosonic quantum 
Hall effect [T^l- This has led to the experimental investigation of a BEC under critical 
rotation , where the rotation frequency equals the trap frequency, and super-critical 
rotation [T?j . 

In this work we investigate the intrinsic reponse of an atomic BEC to rotation 
in a weakly-elliptical trap using the zero-temperature Gross-Pitaevskii equation. This 
is performed over a large range of rotation frequencies, allowing us to probe diverse 
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regimes of condensate dynamics relevant to current BEC experiments. In the regime 
of sub-critical rotation, a quadrupole shape mode dominates the dynamics. Close to 
resonance this mode is unstable and vortex lattice formation occurs, in close agreement 
with experiments. We simulate the role of thermal atoms in the lattice formation, and 
consider the effect of terminating the rotation. At critical rotation, 1] ^ w r , we show 
that large centre-of-mass oscillations dominate the dynamics. 

At ultra-cold temperature T«T C (where T c is the BEC transition temperature) 
and for a large number of condensate atoms, thermal and quantum effects are minimal, 
and an atomic BEC can be approximated by a mean-field 'macroscopic wavefunction' 
ip(r,t). The dynamics of this wavefunction are well-described by the Gross-Pitaevskii 
equation (GPE) Qj)|, 
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where Kxt( r ^) is the external confinement and m is the atomic mass. The nonlinear 
coefficient g = 47rh 2 Na/m governs the atomic interactions, where N is the number of 
atoms and a is the s- wave scattering length. It is useful to consider the hydrodynamic 
interpretation of the wavefunction ip(r,t) = \fn(r, t) exp[z0(r, £)], where n(r, t) is 
the atomic density and 0(r, t) is a phase factor, which defines a fluid velocity via 
v(r,t) = (ft/ra)V0(r,t). 

Vortices in BEC are characterised density node and an azimuthal phase slip, 
quantized to multiples of 2tt (to keep ip single- valued) . The lowest energy state of 
a BEC in a rotating system (above the critical rotation for the onset of vortices) is 
an ordered triangular lattice of singly-charged vortices J20l (vortices with multiple 
'charge' are unstable EB), whereby the system mimics solid-body rotation. However, 
an energy barrier exists for vortices to enter the condensate and so the transformation 
process from a vortex-free state to a vortex lattice is non-trivial. Lattices have been 
formed thermodynamically E2J by condensing a rotating thermal cloud [2HJ . In 
|2H, the precession of a small laser 'stirrer' above a critical velocity generated vortex 
excitations locally [2SJ and ultimately a lattice 10 . Other lattice experiments P3E1IS1 
involve the direct rotation of the condensate in an anisotropic trap, which is the 
case we will study in this paper. In addition we consider a sudden turn on of the 
trap rotation rather than the adiabatic introduction of the rotating anisotropy |BJ . 
Condensate shape oscillations are excited, and through a dynamical instability of the 
dominant quadrupole mode IE] > the condensate fragments [T2] . Provided the 
two-fold rotational symmetry of the system is allowed to break, a turbulent state of 
sound waves (density excitations) and vortices is generated. Subsequently, through 
the interaction of the vortices with the sound field J2Ej , the vortices gradually become 
ordered into a lattice configuration. 

The trap rotation, which occurs in the x-y plane about the z-axis, strongly 
polarizes the system. Although bending and Kelvin wave excitation of vortex lines do 
occur [2jj, the 2D dynamics in the rotating plane dominate, and a 2D description is 
sufficient to simulate lattice formation [H3 E21 E3I • We therefore proceed by solving 
the computationally- advantageous 2D GPE in the rotating x-y plane. 

In the experiments of [2] |EJ the trap is approximately elliptical in the rotating 
plane with the form, 

V ext (x, y, t) = ±mco 2 r [(x 2 + y 2 ) + s(X(t) 2 - Y(t) 2 )] . (2) 
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Here the first term represents the static radial harmonic trap with frequency u r . In 
the second term, e is the ellipticity parameter of the trap (typically £ C 1), and the 
rotating coordinates X(t) and Y(t) are given by, 

X(t) = x cos(fit) - y sin(fit) (3) 
Y(t) = xsm(ttt) +2/cos(fit), (4) 

where Q is the trap rotation frequency. Note that in [3j the trap is strongly 
anharmonic, rather than weakly-elliptical. 

Our units of length and time are the healing length £ = hf^/m/i and (£/c), where 
c = y/i/m is the Bogoliubov speed of sound. The chemical potential /a = nog is 
the unit of energy, where no is the peak atomic density. Unless stated otherwise, 
we employ parameters which match the experiment of Madison et al. 1 j: our trap 
ellipticity is e — 0.025, the chemical potential is \i = lOfkj, the trap frequency is 
uo r = 0.1 (c/£), and the Thomas- Fermi radius of the condensate is Rtf = 14£. For 
such a 87 Rb condensate, our units of time and space correspond to £ ~ 0.2 fim and 
(£/c)~10- 4 s. 

Our initial state is found by propagating the GPE in imaginary time (t — > —it) 
subject to the initial potential V^ xt (x, y : t = 0) while forcing the peak density to unity. 
Subsequently the system is run in real time. We generally characterise the state of the 
system in terms of the condensate energy, which is evaluated according to the GPE 
energy functional, 



E 



I 



dxdy. (5) 



Note that we subtract from equation (5) the energy corresponding to the initial state, 
i.e. energy is measured relative to the ground state. 

In section |21 we investigate the quadrupole mode of a condensate in a rotating 
trap, and discriminate between the stable and unstable regimes. In the latter case, 
we observe vortex lattice formation, and in section [3| we compare our results with the 
lattice experiments. Section BI infers the role of thermal atoms in the lattice formation, 
and the effect of terminating the rotation is analysed in section [5J The condensate 
dynamics under critical and super-critical rotation of the trap are studied in section 
El Finally, section [3 presents the conclusions of this work. 



2. The Quadrupole Mode 

In the rotating frame, and for slow rotations, the stationary solution in an elliptical 
trap is an elliptical condensate EH| > with the size of the ellipticity depending not 
only on the trap ellipticity e but also the rotation frequency ft. Such a state can 
be reached by adiabatically ramping up the rotation frequency from to the chosen 
value, ft 7 j. This technique has been implemented experimentally 8 . In our case, 
however, the rotation is suddenly turned on with frequency Q and maintained at this 
value. The condensate is continuously driven out of equilibrium and never reaches 
a steady elliptical state. Consequently, collective shape excitations are expected to 
occur. 

Classical hydrodynamics predict that, for a non-rotating trapped condensate, the 
shape oscillations satisfy the dispersion relation cj£ = V^W, where £ is the angular 
momentum quantum number [2HJ . The two- fold symmetry of our system implies that 
only even angular momentum modes will be excited, and we expect the lowest order 
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mode to dominate, i.e. the £ = 2 quadrupolar mode with frequency Uq = \^uj r . When 
driven resonantly at half this frequency, i.e. at ft = u r /y/2, this mode is predicted to 
be dynamically unstable [7j, leading to vortex lattice formation. Experimental results 
121 Oj an d simulations of the GPE El lEj have verified this instability, and 

that it occurs over a range of rotation frequencies in the vicinity of coq. 

2.1. Stable quadrupole mode 

We first consider the case of the stable quadrupolar mode. Figure ^a)-(b) shows the 
condensate density and phase throughout a cycle of the quadrupole oscillation for a 
rotation frequency Vt = 0.7uj r . Note that, due to the rotating trap, the two primary 
axes of the quadrupole continuously rotate over time at the rotation frequency Q. 

The condensate density initially has a circular shape. Over time it becomes 
elongated (figure [I] (a) (i)) up to some maximum point (figureC{a)(ii)). The elongation 
then decreases (figure nja)(iii)) and the condensate returns to a circular state (figure 
ma)(iv)). The phase distribution initially, and after one complete cycle, is uniform. As 
the condensate becomes elongated (figure (b)(i)- (hi)), the phase develops a quadrupolar 
structure of the form y) oc XY, where X and Y are the rotating coordinates of 
the ellipse. This corresponds to a flow of fluid towards and away from the trap centre 
in different quadrants of the ellipse, as indicated by the arrows in figure IH[b)(i)- (hi). 
Although the condensate appears to be rotating, the combined effect of the quadrupole 
cycle and the rotation of the quadrupole axes ensures that the fluid remains irrotational 
at all times |Sl[2Hj. When the condensate is highly elongated (figure DJii)), the phase 
structure is dense, i.e. the fluid velocity is high, as required to maintain the irrotational 
flow. 

The condensate energy (figure [He)) oscillates sinusoidally at the quadrupole 
frequency coq between E = 0, representing the circular, ground state condensate, 
and a maximum value E = £ max , corresponding to the fully elongated, excited state. 
This is analogous to a vortex cycling in and out of a condensate driven by a narrow 
optical stirrer jSHj, and to the Rabi oscillations of an atom in a light field. Note that 




Figure 1. Stable quadrupole mode, (a)-(b) Snapshots in density and phase 
during one oscillation of quadrupolar flow (fl = 0.7uj r ), at times t = (i) 185, 
(ii) 370, (iii) 555, and (iv) 740 (£/c). Each box corresponds to the region 
[— 25, 25]£ x [— 25, 25]£. Arrows in (i)-(iii) indicated the direction of fluid flow. 
Density and phase are shown in the ranges (0 — no) and (— 7r, 7r) respectively, with 
light/dark regions corresponding to low/high values. Note that the complex phase 
structure at large radii has no meaning since the density is essentially zero, (c) 
Energy of the condensate during quadrupolar flow. The times of the snapshots in 
(a)-(b) are indicated. 
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Figure 2. The quadrupole mode, (a) Quadrupole frequency uoq as a function 
of rotation frequency Q for the numerical data (circles) and the classical 
hydrodynamic prediction of equation IHl f grev /brown line), (b) Energy amplitude 
£ max of the quadrupole oscillations as a function of rotation frequency. The y- 
values are shown on a logarithmic scale. Inset: the central region of the main plot 
shown with decimal scaling on the y-docis. In all plots the dashed lines indicate 
the extent of the unstable quadrupole region. 



the energy shows no additional behaviour: we do not detect higher order modes (e.g. 
I = 4, 6, ...), and the orientation of the quadrupole (oscillating at Q) does not affect 
the energy. 

2.2. Behaviour of the quadrupole mode with rotation frequency 

The sinusoidal energy evolution in figure [He) is a general feature of the quadrupole 
mode, characterised by the quadrupole frequency uoq and amplitude E max . We now 
examine the variation of these quantities with rotation frequency Q. First, consider 
the quadrupole frequency, shown in figure Efa). In the limit fi — > the quadrupole 
frequency approaches uJq = \/2uj r) which corresponds to the analytic prediction in 
a non-rotating system [221 • F° r a system rotating at frequency Q the Hamiltonian 
features an additional centrifugal term —QL Z 20 , where L z is the ^-component of 
the angular momentum . This adds a shift of — £Q to the dispersion relation. For the 
quadrupole 1 = 2 mode this gives, 

u> Q (Sl) = \cu° Q - 2fi|. (6) 

This equation is plotted in figure ^a) (grey/brown line) and fits the numerical data 
well for rotation frequencies away from the quadrupole resonance at Q = 
However, close to this resonance, the data deviates from the linear trend and never 
reaches zero. This is the unstable region of the quadrupolar mode, which ultimately 
leads to vortex lattice formation . Note that the unstable regime lies in the vicinity 
of, but not necessarily centred around, the analytic prediction of Q = ujq/2. This 
shift is due to the elliptical perturbation on the system. The extent of this unstable 
region and how it compares to experimentals will be discussed in section The 
amplitude of the quadrupole mode, shown in figure EJb), increases towards the region 
of instability. Close to and within this region, the amplitude becomes extremely large. 
The unstable quadrupole region is therefore associated with when the quadrupole 
mode has a frequency close to zero and a large amplitude. 
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3. Comparison with experimental observations 

3.1. Width of the unstable region and number of vortices 

In the experiment of Madison et al. Q], lattice formation was observed in the 
region 0.68 < Q/u r < 0.78. Below this no vortices were observed, while above this 
the condensate became turbulent. Figure |3Ja) shows the final number of vortices 
we observe in the lattice for different rotation frequencies. To be consistent with 
the experimental observations, where density imaging has limited accuracy, we only 
count a vortex if it resides in a region of significant density, which we take to be 
where the density is greater than n « 0.2no- This also avoids the counting of the 
phase singularities at very low density which have negligible energy. The range of 
rotation frequencies over which we observe lattice formation, 0.72 < Q/u r < 0.78, is 
comparable to the experimental observations. Within this region, we observe the final 
number of vortices to increase with rotation frequency up to 12 vortices. 

Both below and above the lattice-forming region we observe stable quadrupolar 
oscillations, with no indication of a turbulent regime. Note that for super-critical 
rotations Q/co r > 1 no more condensate atoms were observed in the experiment. We 
will consider this regime in section 
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Figure 3. Number of vortices in the final lattice versus rotation frequency for (a) 
the system of Madison et al. 1 with trap ellipticity £ = 0.025, and the system 
of Hodby et al. 2 with trap ellipticities of (b) e = 0.041 and (c) e = 0.084. 
Numerical results are show by crosses. In (b)-(c) the experimental data (averaged 
over four observations) from |2] is shown by circles. 



We now consider the experiment of Hodby et al. , and employ a system with 
chemical potential \i « Afiuo r and trap frequency uo r « 0.25(c/£), which has a Thomas- 
Fermi radius of Rtf ~ 6£. Figure EJb)-(c) shows the number of vortices produced as a 
function of Q for two different trap ellipticities of (b) e — 0.041 and (c) e = 0.084. The 
experimental observations of Hodby et al. [2j are shown by circles, while the numerical 
results from the GPE are shown by crosses. We see good agreement between theory 
and experiment, both in the extent of the lattice- forming region and the number of 
vortices produced, and particularly for the higher value of the ellipticity. By comparing 
figures Etb) and (c), we see that the width of the unstable regime and the number of 
vortices produced increase with trap ellipticity. The number of vortices produced here 
is lower than for the parameters of Madison et al. since the condensate is smaller. 
Note that Lundh et al. ^UJ have modelled this experiment using the 2D GPE (albeit 
for different data) and also found good agreement. 
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3.2. Symmetry breaking and the timescale for lattice formation 

In our simulations of the unstable quadrupole regime, numerical noise grows 
exponentially and ultimately breaks the symmetry on a macroscopic level [T21- In 
real systems, noise arises from many sources, e.g. thermal/quantum fluctuations, 
trap imperfections and stray fields, and will have a similar symmetry-breaking effect, 
albeit on a much larger scale. In order to model such 'real-life' effects we impose 
a random 'jitter' on the trap centre, limited to the region [—7,7] x [—7,7]. The 
parameter 7 therefore allows us to control numerically the degree of symmetry- 
breaking. Conveniently, even a very small jitter, e.g. 7 = 0.0001£, is sufficient to 
dominate over the effects of numerical noise 12 . However, the final state of the 
lattice, e.g. the number of vortices, is insensitive to 7. We will employ a jitter of 
7 = 0.1£, which corresponds to around 20 nm. 

We return to using the system parameters of Madison et al. 1 . Figure Ufa) 
shows the evolution of the condensate energy for various rotation frequencies within 
the unstable region. Figure 0{b) shows the corresponding evolution of the elliptical 
deformation of the condensate a defined as, 

<X 2 >-<Y 2 > _ 
a = n <X*> + <Yi> ' (7) 
where < X 2 >= J X 2 \^(x,y = 0)| 2 dX and < Y 2 >= j Y 2 \^(x = 0,y)\ 2 dY. These 
curves have the same qualitative features, which are related to the key stages in the 
evolution [T2| . 

Fragmentation: The quadrupole mode becomes over-excited and atoms are ejected. 
Partial quadrupole oscillations occur (where the condensate does not fully return 
to its circular state), with the energy undergoing partial oscillations. 

Symmetry breaking and turbulence: When the system asymmetry has grown to 
macroscopic levels the quadrupole oscillations cease. Energy is rapidly coupled 
into the condensate up to some maximum value. At this point the condensate is 
a turbulent state of vortices and sound waves. Note it also returns to having an 
approximately circular shape. 

Crystallisation: The total energy stays approximately constant, while vortex-sound 
interactions allow the vortices to settle into a lattice configuration. 

As apparent in figure the length of the fragmentation stage is highly sensitive 
to the trap rotation frequency. In the shortest case (solid black line in figure 0), 
the fragmentation stage persists for a time of around 3000 (£/c) (or 300 ms). On the 
other hand, the timescale for the symmetry-breaking/turbulent stage is insensitive 
to these effects, and typically extends for a time of around 1000 (£/c) (or 100 ms). 
The crystallisation phase in our simulations is also independent of these effects and 
remains a slow process, typically taking 50000 — 10000 (£/c) (or 0.5 — 1 s) to form a 
well-ordered lattice. 

In the experiments El El? the condensate is observed to undergo partial 
quadrupole oscillations (where the condensate does not fully return to a circular state), 
become disordered, and finally form an ordered vortex lattice. This appears to be 
consistent with the stages of lattice formation described above. Indeed, in [HI the 
elliptical deformation of the condensate is measured as a function of time, and bears 
the same qualitative features as in figure ^b). 

We now compare our simulations to the experimentally-observed timescales, 
where possible. In 8 , the partial shape oscillations occur for approximately 250 ms. 
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Figure 4. (a) Evolution of the condensate energy for rotation frequencies 0.73 
(solid grey /brown line), 0.75 (dashed black line) and 0.77 (solid black line), (b) 
Evolution of the elliptical deformation of the condensate, denned by equation J7|, 
for the corresponding cases. A trap jitter of 7 = 0.1£ is employed. Note that the 
time unit corresponds to (£/c) ~ 10 -4 s. 



Subsequently, the experimental shape oscillations cease and the condensate returns to 
an approximately circular shape. This occurs over a time of around 50 ms. These 
timescales are comparable with the quickest evolution in figure 21 (solid lines). 

In the experiments, a lattice is observed after (0.5 — 1) s of trap rotation, 
where in our simulations [E it takes the equivalent of approximately 10 seconds 
for the lattice to become well-ordered. We have shown that the timescale for 
the fragmentation, symmetry-breaking and turbulent phases are consistent with 
experimental observations, and so this time anomaly must lie with the crystallisation 
phase. This discrepancy in the crystallisation rate may be due to thermal dissipation 
(although results in section EI would contradict this) or 3D effects such as Kelvin wave 
excitation and damping [2jj, both of which are not present in our simulations. 

4. The role of 'thermal' atoms 

In the classical field interpretation JSJ, the GPE describes both the mean- field 
condensate and the fluctuating thermal field of random excitations. The thermal 
excitations typically have high energy and momenta in comparison to the condensate. 
As described in > an integral part of the lattice formation process is the irreversible 
generation of a highly-energetic outer cloud, which we will interpret as thermal atoms. 

The numerical grid on which the simulations are performed naturally imposes a 
momentum cutoff into the system. In our square 2D system, with isotropic grid spacing 
the £;-values are restricted to a square region \k x ^ y \ < & max , where /c max = 2tt/S. We 
employ an isotropic grid spacing S = 0.25£, which is sufficiently small that it does not 
restrict directly the dynamics of the condensate. However, it will affect the thermal 
component of the system. High energy atoms will 'evaporate' from the system when 
they exceed this cutoff. For example, consider the total energy and norm during the 
lattice formation process (figure [SJ solid black lines) . The total norm does not visibly 
decrease, yet a very small proportion of thermal atoms are lost. These are the most 
energetic atoms, causing a gradual decrease in the energy at later times (figure |^a), 
solid black line). 

We investigate the role of the classical thermal cloud in the lattice formation 
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Figure 5. Effect of momentum cutoff on lattice formation, (a) Evolution of the 
condensate energy for the unimposed case (grey /brown line), where the allowable 
fc-values are restricted by the numerical grid to a square region [k x ,k y < fcmax], 
and imposed circular regions with radius k cu t = fc max (solid black line), 0.8fc max 
(dashed black line), and 0.6/c m ax (dotted black line). Here fe max = 2tt/5 is the 
maximum allowable /c-value, which is set by the isotropic grid size 5. (b) The norm 
of the condensate for the same cases, (c) Number of vortices in the condensate 
for these four cases. 



process by varying the cutoff in &-space. To impose a cutoff at fc C ut? we Fourier 
transform into &-space, set \ip(k x , k y )\ — for values which satisfy (k x + k 2 ) 1 ^ 2 > fc cu tj 
and then inverse Fourier transform to obtain the modified spatial wavefunction. This 
is repeated regularly throughout the simulation (every 100 time steps). Note that the 
imposed cutoff is circular in &-space, in contrast to the square cutoff induced by the 
grid. 

FigureElalso shows the dynamics for various imposed circular cutoffs: k cut = k max 
(grey/brown solid line), 0.8/c max (grey/brown dashed line), and 0.6/c max (grey/brown 
dotted line). These cases exhibit some common features. Despite the reduction in 
the population of thermal atoms, a vortex lattice forms in all cases. Importantly, the 
lattices contains approximately the same number of vortices (figure Et c ))- The peak 
condensate energy (figure Eta)) is slightly lower than the case of no imposed cutoff 
(black line). In the crystallisation phase the energy decreases exponentially, and at a 
faster rate than the case of no imposed cutoff. Similarly, for the norm (figure Etb)), 
atoms are lost during the lattice formation process, and the rate of loss increases as 
the cutoff is reduced. These features are to be expected since lowering the cutoff 
'evaporates' more energetic atoms from the system. However, the number of vortices 
in the condensate does not visibly decrease, indicating that this energy /particle loss 
has little or no affect on the vortex dynamics in the central condensate. 

What is somewhat surprising is the significant difference in the energy curves 
between the square cutoff (figure Et a )? black solid line) and the circular cutoff with 
the same width (figure Et a )> grey /brown solid line). Clearly a large amount of energy 
must acculumate in the high- A: corners of &-space. 

For the imposed cutoffs, the energy has very similar dynamics, despite the large 
variations in the cutoff. In particular, the exponential decay of the energy during the 
crystallisation stage occurs on similar timescales. Furthermore, the number of vortices 
in the final lattice is insensitive to the cutoff. This suggests that the formation of the 
lattice is independent of the number of thermal atoms, and therefore independent 
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of temperature. Experimental evidence to suggest that the crystallisation stage is 
temperature-independent has been presented in [3 EJ . 

5. Effect of terminating the rotation 

So far we have considered a condensate in a trap which is under continuous rotation. 
Here we study the effect of terminating the rotation after some time. We first consider 
the case of the stable quadrupole mode (described in section (2J). This cycle is driven 
and maintained by the continuous trap rotation. Terminating the rotation freezes the 
energy and elongation of the condensate. The quadrupole mode ceases although the 
elongated condensate stills rotates at the rotation frequency. 




4000 8000 12000 



t &c) 

Figure 6. Effect of terminating the rotation. Evolution of the condensate energy 
for continuous rotation (black line), and rotation which is terminated at times 
t = (i) 1180, (ii) 4000, (iii) 7200, and (iv) 8000 (f/c). The number of vortices 
produced in each case is indicated. 

Figure HO presents results for the unstable quadrupole regime. Here the energy 
for the case of continuous rotation (black line) is compared with cases where the 
rotation is turned off at various times (grey/brown lines). Once the rotation is 
terminated, no more energy can be coupled into the system, and the total energy 
becomes approximately fixed (some energy is lost due to the finite momentum cutoff 
- see section 0J. In all cases shown a vortex lattice is ultimately produced, including 
the case where the rotation is terminated after just one quadrupole oscillation. This 
implies that once the instability has been started, i.e. the fragmentation during the 
first quadrupole oscillation, the full vortex formation process will occur. 

An interesting effect of terminating the rotation is that, because the energy of 
the system becomes fixed, it allows control over the number of vortices in the final 
lattice. For example, under continuous rotation the final state in figure El (black line) 
contains a maximum value of 11 vortices. However, for the cases where the rotation 
is terminated (grey/brown lines) the final state contains 1,2,3 and 4 vortices. In 
principle using this technique it should be possible to generate final lattices containing 
any number of vortices up to the maximum value. 

6. Critical rotation 

The fast rotation of condensates at the critical frequency fi = uo r and above is of 
current interest, with the prospect of ultimately generating exotic phenomena such as 
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vortex melting and a quantum Hall state of bosons [T5|. For a gas of classical non- 
interacting particles under critical rotation the centrifugal force equals the trapping 
force, and the particles are no longer confined 31 j. Super-critical rotation has been 
studied experimentally in a harmonic trap featuring an additional quartic term to 
avoid this singularity [T7|. Rosenbusch et al. [El investigated critical rotation in 
a rotating trap both experimentally and theoretically, and showed that, while the 
non-interacting (non-condensed) gas explodes, repulsive interactions can hold the 
condensate together. This occurs if the trap ellipticity does not exceed a critical value 
(of the order of e ~ 0.1). Then the centre of mass of the condensate is destabilised 
and spirals away from the trap centre. However, if e exceeds the critical value, the 
condensate also undergoes explosive dynamics. Note that the values of ellipticity used 
throughout this paper are well within this critical ellipticity. 

We now consider the response of our system to critical and super-critical trap 
rotation. As shown in figure there is no strange behaviour in the quadrupole mode 
as the critical point is reached and exceeded - the frequency /amplitude continue to 
increase/decrease in a smooth manner. However, we do observe additional dynamics 
under critical rotation, namely the periodic spiralling of the condensate in and out 
of the trap centre. This is in partial agreement with Rosenbusch et al. ^2), who 
conclude that spiralling condensate does not return to the centre of the trap. As with 
lattice formation, these dynamics require the initial two-fold symmetry of the system 
must be broken. Although this will occur slowly due to numerical noise, we impose a 
trap jitter to speed up the dynamics, unless stated otherwise. Figure 0Ja) (solid line) 
shows the evolution of the condensate energy for critical rotation Q = uj r . The energy 
undergoes huge oscillations. This corresponds to the periodic centre-of-mass motion of 
the condensate to large radii in the trap (figure Efb)). Quadrupolar shape oscillations 
of the condensate still occur, but the energy amplitude of the quadrupole oscillations 
(see figure EJb)) is three orders of magnitude smaller than these oscillations. In 
the absence of an imposed trap jitter (dotted line in figure Ufa)), the centre-of-mass 
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Figure 7. Critical rotation, (a) Evolution of the total condensate energy for 
critical rotation = uj r with trap jitter 7 = 0.1 (solid line line), and with no 
trap jitter (dotted line), (b) Radius of the condensate peak (with trap jitter), (c) 
Density and (d) phase plots of the condensate during the centre-of-mass precession 
at times t = (i) 500, (ii) 3500, and (iii) 4500 (£/c). These times are indicated in 

(b). 
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Figure 8. Rotation at fl = 0.99o; r . (a) Condensate energy for a rotation 
frequency Q = 0.99uu r (with trap jitter 7 = 0.1£). Insets show the path of the 
condensate peak during the times (i) 3000 - 6000, (ii) 9000 - 12000, and (iii) 
18000 — 23000 (£/c). (b) Initial energy amplitude £ max of the centre-of-mass 
motion. 



dynamics begin after the system has evolved for some time at low energies. 

FigureUb)-(c) shows the density and phase during the centre-of-mass oscillations. 
The condensate begins with approximately circular shape and uniform phase (snapshot 

(i) in figure |7fc) -(d)). It slowly drifts away from the trap centre (figure Efc), snapshot 

(ii) ) while precessing quickly at the rotation frequency. A striking quasi-linear phase 
structure develops (figure [7fd), snapshot (ii)), which represents the bulk motion of the 
condensate in the tangential direction. At a later time (figure [TJc)-(d), snapshot (hi)), 
the condensate is further from the trap centre, and here the linear phase structure is 
more concentrated. Periodically the condensate returns back to the trap centre. 

We observe that this centre-of-mass motion occurs in a range of frequencies around 
the critical frequency. However, the centre-of-mass precession is not always of constant 
amplitude and frequency. Figure Efa) shows the condensate dynamics for a rotation 
frequency of Q = 0.99o; r . Here the energy oscillates with varying frequency and 
amplitude. As shown in the insets of figure Efa), this is due to the condensate moving 
in an ellipse, rather than a circle, with the ellipticity and axes of the ellipse clearly 
changing with time. 

We characterise the size of the centre-of-mass motion in terms of the amplitude 
of the initial energy oscillations. Figure Efb) shows this energy amplitude E max for 
rotation frequencies around the critical point. We clearly see that the centre-of-mass 
motion occurs in the narrow frequency range 0.985 < Q/uj r < 1.015. In this range 
is predicted to be y/\ — £ < Vt/uo r < y/1 + i, giving 0.987 < Vt/uo r < 1.012, which is 
in good agreement with our observations. The irregular centre-of-mass motion occurs 
in the lower part of this range, for Vt/uo r < 0.9975, where the condensate is excited to 
the highest energies. 

Above this critical region, the centre-of-mass oscillations cease and the dominant 
excitation in the condensate returns to being the quadrupole mode, albeit with very 
low amplitude (see figure EJ). 

We have also considered the critical rotation of a condensate which contains 
vortices, both by forming a vortex lattice at ft ~ 0.7w r and suddenly switching the 
rotation to ft = uo ri and by slowly ramping the rotation, through the unstable regime, 
up to Vt = uo r . In both cases, we observe the condensate to undergo centre-of-mass 
oscillations in the trap. However, here the oscillating condensate typically has a diffuse 
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and disordered appearance, and contains several vortices. 

7. Conclusions 

In this paper, we have investigated the dynamics of an atomic Bose-Einstein 
condensate in a weakly-elliptical rotating trap using the two-dimensional Gross- 
Pitaevskii equation. This was performed over a wide range of rotation frequencies 
relevant to current experiments. The quadrupole mode of the condensate plays a key 
role in the dynamics, with the fate of the system depending crucially on the rotation 
frequency. We resolve the stable and unstable regimes of the quadrupole oscillations. 
The unstable regime occurs around the resonant frequency of the quadrupole mode, 
cJq = (j r /V% an d rotation within this range ultimately leads to vortex lattice 
formation. Indeed, over the large range of rotations investigated here, we only observe 
vortex nucleation to occur within this range, i.e. vortex nucleation in elliptical traps 
only occurs through the dynamical instability of the quadrupole mode. We compare 
with lattice experiments and find good qualitative, and sometimes quantitative, 
agreement, despite the limitations of the model we employ (two-dimensional and 
zero-temperature). This reinforces the idea that lattice formation is at heart a two- 
dimensional effect. Furthermore, we find that the process is insensitive to the cutoff of 
'thermal' atoms in the system, suggesting that lattice formation is a zero-temperature 
phenomenon. 

By terminating the trap rotation during the lattice formation process we show 
that it is possible to control the number of vortices in the final lattice. Finally, 
we considered the response of the system to critical rotation at the trap frequency, 
and observe large centre-of-mass oscillations of the condensate, without the explosive 
dynamics expected for a classical particle. 
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